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FEEDBACK  CONTROL  OF  LINEAR  TURBULENCE  USING 
ELECTROMAGNETIC  MICROTELES 


1.  INTRODUCTION 

Control  of  turbulence  leads  to  reduction  in  the  viscous  drag  and  suppression  of  turbulence- 
induced  noise.  Via  direct  numerical  simulation  of  a  low  Reynolds  number  channel  flow  over  a 
riblet  surface,  Crawford  and  Kamiadakis1  have  shown  that  drag  reduction  is  uniquely  related  to 
the  suppression  of  the  surface-normal  component  of  turbulence  near  the  wall.  Bushnell,2 
Bandyopadhyay,3  and  Gad-el-Hak4  describe  research  done  in  this  area.  Several  experimental  and 
computational  studies  have  also  focused  on  the  delay  of  boundary-layer  transition  through  wave 
suppression  by  introducing  waves  of  appropriate  amplitudes  and  phases.5,6,7  In  recent  studies, 
Joshi8  and  Joshi  et  al.9  have  taken  system-theory  approaches  to  the  channel  flow  control  by 
suction  and  injection  at  the  wall.  In  a  system-theory  approach,  the  model  is  initially  represented  by 
a  set  of  first-order  differential  equations,  and  then  modem  system  and  control  theory  is  applied  to 
the  model  design  and  to  system  analyses. 

A  Lorentz  force  field  is  produced  when  electric  and  magnetic  fields  are  applied  in  an 
electrically  conducting  medium.  The  potential  application  of  Lorentz  forces  for  drag  reduction 
and  turbulence  control  has  been  investigated  by  Henoch  and  Stace10  and  Nosenchuck  and 
Brown.11  These  approaches  to  electromagnetic  turbulence  control,  however,  do  not  use  feedback 
algorithms.  Nosenchuck  and  Brown11  have  pulsed  a  surface-normal  Lorentz  force  in  an  attempt  to 
minimize  the  outward  component  of  the  surface-normal  perturbations.  On  the  other  hand,  Meng12 
has  proposed  a  sensor-based  active  feedback  control  scheme  and  has  modeled  the  turbulence 
production  process  as  the  appearance  of  a  sequence  of  organized  motions  given  by  the 
probabilistic  distributions  of  their  scales.  Meng  then  proposed  a  Markov  chain  control  scheme  that 
uses  sensors  and  attempts  to  control  the  Lorentz  forces  over  all  phases  of  the  turbulence 
regeneration  process.  At  the  Naval  Undersea  Warfare  Center  (NUWC)  Division,  Newport,  RI, 
Bandyopadhyay  has  designed  two-dimensional  thin  arrays  of  electromagnetic  tiles  with  spatial 
scales  on  the  order  of  millimeters  that  are  suitable  for  high  Reynolds  number  turbulence  control.10 
The  electromagnetic  forces  produced  by  these  microtiles  are  confined  to  within  1  millimeter  from 
the  wall  and  are  promising  for  the  boundary-layer  turbulence  control  at  high  Reynolds  numbers. 14 

Low  Reynolds  number  laminar  flow  experiments  have  indicated  that  the  NUWC  microtiles 
produce  “pillows”  of  vorticity  that  scale  with  their  size.  It  is  believed  that  the  microtiles 
redistribute  the  prevailing  vorticity  in  the  shear  flow  but  do  not  introduce  vorticity  into  the  flow. 
The  observations  of  the  pillows  of  vorticity  have  been  later  confirmed  by  direct  numerical 
simulation.15  The  microtiles  can  be  described  as  hairpin  vortex-generators.  Further  numerical 
simulations  of  a  turbulent  channel  flow  over  an  array  of  two  microtiles  at  low  Reynolds  numbers 
and  subjected  to  a  steady  power  supply  indicate  that  an  8-percent  suppression  of  turbulence  can 
be  achieved  by  the  low-power  microtiles.  16In  this  simulation,  the  channel  flow  Reynolds  number 
Re  was  50,470.8,  based  on  channel  half- width  and  maximum  centerline  velocity.  The  magnets  and 
electrodes  in  each  microtile  were  arranged  in  a  checkerboard  pattern,  where  one  ideal  Lorentz 
force  was  vertical  (in  reality,  one  force  is  three-dimensional  near  one  comer).  Within  a  microtile, 
the  open  area  between  the  magnets  and  electrodes  measured  4  mm  x  1  mm  in  the  streamwise  and 
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cross-stream  directions,  respectively.  The  applied  maximum  magnetic  field  strength  Bo  was  0.6  T 
and  the  voltage  £  was  6.0  V.  An  array  of  2-x-2  microtiles  was  considered,  and  the  turbulence 
suppression  reported  above  was  attained  at  12  wall  units  from  1  surface. 

The  present  study  focuses  on  the  control  of  a  two-dimensional  turbulent  boundary  layer  on  a 
flat  plate  and  explores  the  feasibility  of  a  closed-loop  control.  The  system- theory  approach  taken 
here  is  similar  to  that  of  Joshi,8’9  however  the  Lorentz  force  field  produced  by  microtiles  is  the 
control  mechanism.  Moreover,  a  flat-plate,  boundary-layer  flow  is,  because  of  the  nature  of 
boundary  conditions,  relatively  complicated  compared  with  the  channel  flow.  Although  turbulence 
control  for  any  finite  perturbation  is  the  final  goal,  this  study  is  limited  to  control  of  small 
perturbations  in  the  base  flow.  The  present  approach  to  control  via  Navier-Stokes  equations  is 
rooted  in  the  structural  modeling  of  a  turbulent  boundary  layer  as  advanced  by  Perry  and  Chong17 
and  by  Bandyopadhyay  and  Balasubramanian.18, 19  Structural  modeling  of  a  turbulent  boundary 
layer  is  better  suited  for  feedback  control  than  is  statistical  turbulence  modeling.  In  conventional 
turbulence  modeling,  a  turbulent  flow  is  described  by  the  Reynolds  decomposition,  where  it  is  said 
to  be  composed  of  a  time-mean  flow  and  random  perturbations.  The  time-mean  flow  is  described 
statistically  via  closure  models.  However,  in  this  report,  along  the  lines  of  Bandyopadhyay  and 
Balasubramanian,18’ 19  a  turbulent  boundary  layer  is  described  to  be  composed  of  a  laminar-like 
base  flow  and  random  perturbations.  Because  feedback  control  requires  deterministic  signals,  the 
perturbations  are  described  via  Fourier  series  and  not  statistically.  To  date  there  has  been  no  study 
on  feedback  electromagnetic  turbulence  or  any  control  that  could  probably  achieve  a  higher  level 
of  turbulence  suppression  than  an  open-loop  control  can  achieve. 

The  problem  of  turbulence  control  is  complex  because  the  governing  equations  of  motion 
are  nonlinear  and  infinitely  dimensional  and  control  theory  for  infinite-dimensional  systems  is  not 
well  developed.  In  this  report,  a  finite-dimensional  approximate  model  is  derived  based  on  the 
Galerkin  procedure  and  Chebyshev  polynomials.8,9,20’21,22  Based  on  this  model,  linear  feedback 
control  laws  are  derived  using  optimal  control  techniques  and  also  by  using  a  simple  feedback  of 
wall-shear  stress  for  the  stabilization  of  the  perturbed  flow  to  the  base  flow.  It  is  shown  that  the 
unstable  modes  can  be  controlled  by  longitudinal  as  well  as  surface-normal  forces.  However,  for 
lower  wave  numbers,  longitudinal  forces  are  more  effective  and  require  smaller  electrode  voltages 
for  control  compared  with  surface-normal  forces.  It  is  shown  that  a  judicious  choice  of  spatial 
distribution  of  Lorentz  forces  is  important  for  designing  stabilizing  controllers,  and  this  requires 
proper  activation  of  the  electrodes.  Numerical  results  are  obtained  to  show  that,  in  the  closed- 
loop  system,  asymptotic  stabilization  of  the  flow  is  accomplished  using  moderate  input  voltage  for 
small  perturbations. 

This  report  is  organized  as  follows.  Section  2  presents  the  mathematical  model  and 
addresses  the  control  problem.  A  linear,  finite-dimensional  model  is  obtained  in  section  3.  Section 
4  presents  control  laws,  and  simulation  results  are  given  in  section  5. 
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2.  MATHEMATICAL  MODEL  AND  CONTROL  PROBLEM 


A  two-dimensional  flow  on  a  flat  plate  of  length  l  is  considered.  The  x-axis  is  in  the  flow 
direction,  and  the  y-axis  is  normal  to  the  plate.  The  free  stream  velocity  is  t/«.  Reference  values  h 
=  (8m  /2),  Um,  and  ( h/Ux )  are  chosen  for  the  normalization  of  length,  velocity,  and  time,  where  8m 
is  the  boundary-layer  thickness  (location  where  99.5  percent  of  the  ffeestream  speed  is  reached). 
Suppose  that  r  independent  voltages  are  applied  to  the  microtiles  for  control.  It  has  been  found 
that  Lorentz  forces  are  proportional  to  the  voltages  applied  to  the  electrodes  because  permanent 
magnets  are  used. 

The  nondimensionalized  Navier-Stokes  equations  of  motion  for  the  two-dimensional  flow 
are  given  by 

( dujdt )  +  uipujdx )  +  vidufdy)  -  -{dpjdx)  +  Re~]V2u  +  k0F} (x, y, E ) 

{p v/ dt )  +  u(dv / dx)  +  v(<?V I cfy)  =  -[p p / dy )  +  Re~'V2v  +  k0F2(x,y,E) ,  ( 1 ) 

and  the  continuity  equation  is 
[du/dx]  +  [dv/dy]  =  0, 

where  k0  =  hj{  pU2 ) ;  u  and  v  are  the  longitudinal  and  lateral  velocities;  p  is  the  pressure  field, 
p  is  the  density  of  the  fluid;  Re  is  the  Reynolds  number,  50,470.8;  F\  and  F2  are  the  longitudinal 
and  lateral  Lorentz  forces;  and  EeRr  is  the  vector  of  applied  voltages  to  electrodes.  At  the  wall,  y 
is  -1,  and 

y  =  1  at  the  outer  edge  of  the  boundary  layer  at  x  =  l\  =  dh. 

For  the  basic  flow,  approximate  solutions  of  the  boundary  layer  are  given  as 

u  =  U(y)  =  Uccg(  r|), 
ti  =  (y  +  l)/i/5(x), 

where  g(r|)  is  assumed  to  be  one  of  the  following  forms:25 

si1?)  =  sin(^  77/2) 
or 

giw)  =  2tj  -  2tj3  +  Tj4 . 

Although  the  boundary  layer  on  a  flat  plate  is  a  function  of  x  and  j',  to  simplify  model  derivation,  a 
base  flow  of  the  form 

u  =  U(y),v  =  0,p  =  P(x) 

is  assumed.  This  is  obtained  by  replacing  5(x)  with  8m ,  a  constant  in  the  definition  of  r\.  For  the 
perturbed  motion. 
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u  =  U{y)  +  ux(x,yj),v  =  v,(x,y,f),/?  =  P  +  p(x,y,t), 


(2) 


where  u\,  vh  and  p  represent  small  perturbations  in  the  base  flow.  It  will  be  seen  later  that,  when 
the  base  flow  is  chosen  to  also  be  a  function  of  x  for  linearization,  a  model  is  obtained  in  which  all 
the  modes  are  coupled — creating  complexity  in  controller  design.  Restricting  the  region  from  the 
leading  edge  of  the  plate  causes  the  base  flow  to  be  weakly  dependent  on  x.23 

By  substituting  the  values  in  equation  (2)  in  equation  (1),  neglecting  higher  order  terms, 
eliminating  p\  (by  subtracting  partial  derivative  of  the  second  equation  with  respect  to  x  from  the 
partial  derivative  of  the  first  equation  with  respect  toy  in  equation  (1)),  and  using  the  continuity 
equation,  the  linear  dynamics  are  obtained  by 

[d2ujdtdy)-{d2vjdtdx)  +  Gr(y)((<?2w1/<?x<?'y)  -[d2vjdx2^  +  £/t2)Vj 
=  Re~x{[dzuJ  dx2dy)  +{d2ujdy2)-{d2vjdxdy2)-{d:,vjdxi^) 
+K((eFJffy)-(SFJSx)),  (3) 

where  U'2!  =  d2u(y)/dy2 ,  and  the  arguments  of  various  functions  are  suppressed  for  simplicity. 
The  continuity  equation  for  the  perturbations  is 

(dujdx) +(dvjdy)  =  0.  (4) 

Theoretically,  the  boundary  conditions  are  such  that  u\  and  Vi  tend  toward  zero  as  y  tends  toward 
infinity.  However,  the  effect  of  viscosity  is  negligible  beyond  the  outer  edge  of  the  boundary  layer. 
Therefore,  in  order  to  obtain  a  simplified  model,  the  boundary  conditions  are  assumed  to  be 

K1(-l)  =  V,(-l)=0, 

w,(l)  =  Vj(l)  =  0.  (5) 

The  problem  is  to  design  a  feedback  control  system  so  that,  in  the  closed-loop  system, 
perturbations  U\  (x,  y,  0)  and  Vi  (x,  y,  0)  at  t  -  0  asymptotically  converge  to  zero  as  t  tends  toward 
infinity.  In  this  case,  small  but  arbitrary  perturbations  in  velocities  u  and  v  decay  to  zero,  and  the 
perturbed  flow  tends  toward  the  base  flow  asymptotically. 
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3.  A  FINITE  DIMENSIONAL  MODEL 


In  this  section,  a  finite  dimensional  model  is  obtained  using  the  Galerkin's  method  with  the 
Chebyshev  polynomials  as  basis.8, 20-22  The  choice  of  these  Chebyshev  polynomials  has  given 
useful  results  in  the  channel  flow  control  problem  of  Joshi.8,9  Some  other  independent  set  of 
functions  can  be  used  for  derivation.  However,  the  convergence  property  of  the  derived  model  to 
the  actual  system  depends  on  the  choice  of  basis. 

To  this  end,  a  perturbed  stream  function  y/x  (x,y,  t )  is  chosen  for  the  system  (equation  (3)) 
such  that 

",  =dVi{x,y,t)/dy, 

v,  =-d¥\{x,y,t)/dx  (6) 

satisfy  equation  (4).  Substituting  the  values  in  equation  (6)  in  equation  (3)  provides  a  single, 
partial  differential  equation  given  by 

(<?3  y/Jdtdx2)  +  [d V,  /  dtdy1 )  +  u{y)((d2  \j/x  I  dxdy1 )  +  (d3  y/x  /^xJ)) 

-£/(2)  (dy/x  l&x)  =  Re-1  V2  (v  V, )  +  £0  ((dF:  /  dy)  -  (dF2  /  dx)j ,  (7) 


where 


SJ2^/dx\d2/dy2). 

Motivated  by  the  choice  of  periodic  function  of  x  in  the  derivation  of  the  Orr-Somerfield 
equation22  and  due  to  the  periodic  fluctuations  of  the  velocities  in  the  longitudinal  direction,  a 
hypothesis  can  be  formulated  for  a  series  solution  for  the  stream  function  of  the  form 

co  so 

=  ZZ(«™Wcos(7?a0x)  +  bnm{t)^r{na,x))Lm(y),  (8) 

n=0  m=  0 

where  Lm(y),  (m  -  0,  1,. .  .,oo)  are  the  set  of  basis  functions  formed  by  a  suitable  combination  of 
Chebyshev  polynomials  that  satisfy  the  boundary  conditions  (equation  (5)).  In  view  of  equations 
(5),  (6),  and  (8),  it  easily  follows  that  functions  Lm(y )  must  satisfy 


£„(-!)  =  i.d)  =  0. 

a.„(-\)l3y=cLmb)l8y=Q-  (9) 

Here  the  fundamental  wave  number  is  ao  =  (2rc/^,);  «a0  is  the  wave  number  of  the  nth  harmonic 
term  in  the  series  solution,  b0m  =  0;  and  anm  and  bnm  are  the  time-dependent,  generalized 
coordinates  (amplitudes  of  modes). 

Substituting  equation  (8)  in  equation  (7)  and  truncating  the  series  to  obtain  a  finite 
dimensional  model  gives 
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nm  cos[na0x)  +  b„m  sin(«a0x)  ji.L{^(y)  -  ( na0x )2 Lm{y)\ 

«=0m=0V  ' 

=  sin(«a0x)  +Z>„m  cos(«ar0x))»a0(((«a0): V  +  J/(2))zm  W  -  t/Z,(m2)  (y)) 

+^"’  (a.  cos(«g:0x) + bnm  sin(«a0x)) 

((wao  )4 Lm{y)  + -  l{na0xf  +  k0[dFx  / dy~dF2  / dx),  (10) 

where  an  overdot  denotes  differentiation  with  respect  to  time,  and  =  d'Lm{y)/dy' .  In  this 

derivation,  it  is  assumed  that  (M+  1)  basis  functions  LJy)  and  N  sinusoidal  functions  give  a 
reasonable  approximate  model.  If  the  design  with  chosen  values  for  M  and  N  is  unsatisfactory  and 
the  residual  modes  in  the  closed-loop  system  are  not  small,  then  the  controller  needs  to  be 
redesigned  using  a  model  with  larger  values  for  M  and  N. 

It  has  been  shown  in  reference  13  that  for  each  volt  applied  to  the  electrodes,  the  Lorentz 
force  can  be  accurately  expressed  as  ki  exp (-q0  d),  where  d  is  the  lateral  distance  from  the  wall 
and  it  is  negligible  for  d  >  2  mm.  To  this  end,  spatially  distributed  Lorentz  forces  are  expressed  in 
a  Fourier  series  of  the  form  (/  =  1,  2) 

N 

F,=kZ  (Fnn  cos (na0x)  +  Fnil  sin(wa0x))  exp  (qy)E ,  (11) 

n=0 

where  y  €  [- 1 ,  ym  ],  k  =  k0  kx  exp(g),  q  =  -q0  h,  and  Fmk  is  an  r-row  vector.  It  is  assumed  that  the 
Lorentz  force  is  negligible  forjy  >ym.  The  values  of  qo  and  k\  are  obtained  from  the  Lorentz  force 
plot  given  in  reference  13. 

Because  the  trigonometric  functions  form  an  orthogonal  set,  substituting  equation  (1 1)  in 
equation  (10),  multiplying  equation  (10)  by  cos (nceox)  (sin(«<%x)),  and  integrating  over  [0,4]  give 
the  differential  equations  for  anm  (bnm)  given  by 

lianm[l}l\y)-{na,f  Lm(y)^=  Yubnmnaj^(na0)2  U  +  U{2)]Lm(y) -CZZg’M) 

m— 0  Tn^Q 

Lm(y)  +  L^(y)  -  2 (na0)2  L(2)  (>>)) 
+%4i  i  -  ncc0Fn22 )  exp (qy)E , 

£(*-»( 4?  W  - (na0 ) : 2  Lm W)]  =  2 -anmn<*o(((nao Yu  +  U(2])Lm(y)-  UL(2) (7)) 

+Re~'bnm[(ncxX Lm(y)  +  L(*](y)  -  2 (na0)2  L{2) (7)) 

+k{$Fn\2  +nCC0Fn2\)QXV{Ciy)F  ■  02) 

The  inner  product  of  Ls(y)  and  LJy)  is  defined  as22 

Lm  =  (4M, 4W)  =  \+'Ls{y)Lm(y){i-y2y  5dy .  (13) 
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Define  matrices  as  (5  =  0,1,...,M;  m  ~  0,1,...,  M) 


A.  =  ((a.)  J  =  (a(4,2)  -(»« of  a)), 

A,  =  ((A,)  J  =  (a((««  „)4  A +A4)  -  2(  V  0)2  A2’)). 

A2  =((Aj)_) +((»«  „)2c/+c/2)a)} 

4  =  ((4 .  exp(®/)),  •  •  • ,  (lm ,  exp(gy)))r , 

=  kLy{qF ,11  ~na0Fn22) , 

Bn2  =  kLy(qFnU  +na0Fn2l) ,  (14) 


where  (D„,)sm  denotes  the  s — mth  element  of  matrix  Dm  ,  T  denotes  transposition,  and  in  the 
computation  of  Ly  integration  is  performed  only  over  [-1  ,ym]  where  the  Lorentz  forces  are 
significant. 

By  using  the  definitions  of  matrices  in  equation  (14),  equation  (12)  is  written  in  a  matrix 
notation  as 


4o a"  ~  Dnlan  +  Dn2bn  +  BnXE, 

4o  bn  =  Dn2an  +  DJn  +  BnlE ,  (15) 

where  an  -  (an0, ■  ■  ■  ,£W)r  and  b„  =  (b b „iFf .  The  state  subvector  x„  =  (a„T,  b„T)T  is  associated 
with  the  wave  number  nao  and  it  is  described  by 

xn  =  Anxn  +  BnE ,  for  n  *  0,  (16) 

where 

J  d;'0d„  a; A.' 

-  I-A0A2  AoaJ’ 

A=((A»4,)r.(AX2)rf. 

For  n  =  0, 

x0  =  Aox0+B0E,  (17) 


where  x0=a0=  (a00, . . . ,am)T,  A0  =  £>00'Z)0 , ,  and  B0  =  D0lQBm .  The  vector 

*a=(*o,— ,4)  =  (4,4,4;- 
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is  called  the  state  vector  of  the  system  (equations  (16)  and  (17)),  because  xa  completely  describes 
the  perturbed  velocity  vector  in  the  boundary  layer.  For  a  given  initial  condition  xa(0)  and  the 
electrode  voltage  E(t),  equations  (16)  and  (17)  can  be  solved  for: 

x,(t),i  =  0,1,...,  TV,  t>  0. 

These  time-dependent  coefficients  are  substituted  in  equation  (8)  to  obtain  the  perturbed  stream 
function  and  the  velocity  vector  field. 

It  is  interesting  to  note  that  in  view  of  equation  (16),  the  dynamics  associated  with  various 
wave  numbers  are  decoupled  for  E  =  0.  The  decoupling  of  modes  is  not  possible  if  x-dependent 
base  flow  is  taken.  It  will  be  seen  that  decoupling  considerably  simplifies  control  system  design. 

Now  the  problem  of  flow  control  is  solved  by  deriving  feedback  control  laws  of  the  form 

E  =  -Kxc, 

such  that,  in  the  closed-loop  system,  xa(t)  asymptotically  tends  toward  zero,  where  xc  is  a  chosen 
subvector  of  the  state  vector  xa . 
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4.  CONTROL  SYSTEMS 


Using  equations  (16)  and  (17),  the  complete  system  can  be  written  in  a  state  variable  form: 

=  Aaxa  +  BaE ,  (18) 

where  Aa  =  diag (A0,  A\,...,AN)  and  the  input  influence  matrix  is  Ba  =  [btq  , B( , . . . ,BTN ) 7 .  To 
determine  the  fluctuations  of  the  wall-shear  stress  at  a  point  x*,  because  it  is  proportional  to 
du\ldy ,  consider  an  output  variable ya  =  dui/dy(x*,  y=-\,  f)  written  as 

ya=caxa,  (19) 

where  the  output  matrix  Ca  is 

Qj  —  (CaO>-"iCaw)> 

4  =  (4?>(-l),.. 

Ca0  ~  Le, 

and  for  n  *  0, 

can  =(cos(mz0x*)4,sin(«a0x*)4). 

This  output  can  be  used  for  feedback  for  stabilization  or  to  directly  control  ya. 

The  concept  of  controllability  plays  an  important  role  in  system  theory  .24  The  system 
(equation  (18))  is  said  to  be  controllable  if,  given  any  two  arbitrary  points  xa  and  x\  in  the  state 
space,  there  exists  a  control  E  that  can  steer  the  system  (equation  (18))  from  x'a  to  x]  in  a  finite 
time.  The  system  is  controllable  if,  and  only  if,  the  controllability  matrix  CO  defined  as 

co=(bab,...,as-]b), 

has  rank  5  =(2/V+l)(M+l),  which  is  the  dimension  of  state  vector  xa.  The  spatial  distribution  of  the 
Lorentz  forces  must  be  chosen  to  satisfy  the  controllability  conditions.  For  the  controllable 
system,  a  linear  feedback  law  to  stabilize  the  system  can  be  chosen. 

For  clarity,  it  is  assumed  that  the  electrode  voltage  E(t)  is  a  scalar  function.  The  design  for 
the  multi-input  system  can  be  similarly  done.  It  is  interesting  to  note  that  in  the  model  (equation 
(18)),  not  all  modes  are  unstable.  For  the  chosen  parameters  identified  in  section  5,  it  is  found  that 
only  matrix  A\  associated  with  the  fundamental  wave  number  ao  is  unstable,  the  remaining 
matrices  A,  are  stable.  Because  the  electrode  voltage  is  limited,  it  is  a  good  idea  to  design  a 
control  system  to  stabilize  only  the  dynamics  of  unstable  wave  numbers  and  to  leave  the  responses 
of  other  modes  unaltered.  The  dynamics  of  the  unforced  system  are  decoupled,  therefore,  one 
way  to  preserve  the  dynamics  of  the  stable  modes  is  to  choose  the  spatial  distribution  of  the 
Lorentz  forces  such  that  the  control  influence  matrix  Bk  of  the  stable  modes  are  null;  that  is, 
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B0  =  0,  Bn  =  0,  n*  1. 


(20) 


For  this  choice,  when  the  feedback  loop  is  closed  to  stabilize  the  unstable  mode  using  the 
measured  signal^,  control  input  will  not  excite  the  remaining  modes. 

The  condition  in  equation  (20)  can  be  satisfied  by  choosing  F\  and  F2  of  the  form  (z  =  1,  2) 
F,  =  *,[/=;„  cos  a0x  +  FU2  sina0x]exp(©/)£ ,  (21) 

.  ,y 

where  (FJ2,  +F22)  2  =  1 .  The  sinusoidal  distribution  of  Lorentz  forces  requires  appropriate 
activation  of  the  electrodes.  If  other  harmonic  terms  are  present  in  equation  (21),  the  dynamics  of 
other  wave  numbers  will  be  excited  when  the  measured  signal  ya  is  fed  back. 

Because  either  the  longitudinal  or  the  surface-normal  force  for  control  can  be  chosen,  a 
question  arises  concerning  the  effectiveness  of  Lorentz  forces  for  control.  In  view  of  the  input 
influence  matrices  Bn\  and  Bn2  in  equation  (14),  if 


na0>-q. 


the  input  influence  matrix  B„  is  larger  when  surface-normal  force  is  used.  However,  for  smaller 
wave  numbers  not  satisfying  the  above  inequality,  the  longitudinal  force  is  preferable  for 
stabilization  because  smaller  electrode  voltage  will  be  needed  for  control.  In  section  5.2  a  simple 
feedback  law  using  wall-shear  stress  (or  du\/dy(x*,  y  =  -1,  t)  for  stabilization  is  simulated. 

Several  control  techniques  for  linear  systems  are  available  for  the  control  system  design  and 
can  be  applied  to  obtain  different  response  characteristics  and  to  meet  design  specifications.  For 
the  choice  of  Lorentz  force  as  given  in  equation  (21),  the  optimal  control  technique  can  be  used. 
Consider  the  dynamics  associated  with  the  fundamental  wave  number 

Xi  =  AjX}  +  B]E,  (22) 


which  is  unstable.  Here  x,  =  )  is  the  state  vector  associated  with  the  wave  number  ao. 

The  optimal  control  is  obtained  by  minimizing  a  quadratic  performance  index  of  the  form 


J  = 


Ox  j  +  AE2 


0 


(23) 


where  0  is  a  positive  definite  symmetric  matrix,  and  X  >  0.  The  matrices  0  and  A  are  selected  to 
weigh  the  relative  importance  of  performance  measures  caused  by  the  state  vector  perturbations 
x\  and  the  electrode  voltage  E(t).  The  optimal  control  law  is 

E  =  -ZlB(Rx,=-Kxl,  (24) 

where  the  positive  definite  symmetric  matrix  R  is  the  solution  of  the  Riccati  equation24 

RAX  +  AfR  -  RB,A~XBTXR  +  Q  =  0, 
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and  K  is  a  row  vector  of  dimension  2(M+1). 

Substituting  linear  control  law  (equation  (24))  in  equation  (22)  gives  the  closed-loop 
system: 

x,  =(4  ~B,k)x„ 

for  which  the  matrix  ( A ,  -  B}K)  is  Hurwitz  (i.e.,  all  the  eigenvalues  of  (t,  -  BXK)  have  negative 
real  parts).  Thus,  in  the  closed-loop  system,  x,  {t)  asymptotically  tends  toward  zero,  which 
implies  that  the  perturbations  ux  and  vi  decay  to  zero  as  well.  The  response  characteristics  for  the 
state  vector  xx  depend  on  the  eigenvalues  of  the  matrix  (  A]  -  BXK) . 

Substituting  Oyt ,  O  >  0  in  place  of  xf  Ox]  in  the  performance  index  can  be  done  for 
minimization,  where 

Ti  =(cos(a0^*)Zc,sin(a0x*)Zc)x1Ac]x1  (25) 

for  directly  influencing  du\/dy  at  x  =  x*.  By  a  proper  selection  of  the  weighting  matrix,  0,  and  the 
scalar,  X,  a  tradeoff  between  the  speed  of  convergence  of  the  perturbed  state  to  the  equilibrium 
state  and  the  magnitude  of  the  electrode  voltage,  which  is  required  for  control,  can  be  obtained. 
Figure  1  shows  the  complete  closed-loop  system,  including  the  optimal  controller  (solid  line)  and 
the  simple  controller  using  wall-shear  stress  feedback  (dotted  line).  In  the  feedback  loop  using 
wall-shear  stress  feedback,  the  gain  K  is  scalar. 


LINEAR  BOUNDARY 


NOTES: 

SOLID  LINE  -  OPTIMAL  FEEDBACK  LOOP 

DOTTED  LINE  -  WALL-SHEAR  STRESS  PERTURBATION  FEEDBACK  LOOP 


Figure  1.  Closed-Loop  Lorentz  Force  Control  System 
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5.  SIMULATION  RESULTS 


In  this  section  the  results  of  simulation  using  MATLAB®  and  Simulink™  software  are 
presented.  To  derive  the  model,  it  is  assumed  that  the  normalized  velocity  is 
u(y)  =  sin(;r(_y  + 1)/4).  The  parameters  chosen  for  high  Reynolds  number  flow  are  Um  =  17(m/s), 
h  =  (6.77/2)mm;  kinematic  viscosity  of  water,  v  =  1 1 .4£-7m2/s,  v  =  pp,  i\  =  2tc;  and  density  of 
water,  p  =  1000(kg/m3).  For  simulation,  the  values  kx  =  45(N/mJ)  and  qo  =  (l 000/1  -  3)(nf')  are 

used.  These  parameters  are  obtained  from  reference  13.  The  fundamental  wave  number  is  a0  =  1. 
For  a  choice  ofM=  2  in  the  series  expansion  (equation  (10)),  a  sixth-order  state  variable  model 
for  each  wave  number  is  obtained.  The  three  elements  of  the  basis  taken  from  reference  8  are 

A.W=3-4?W+?W, 

L,(y)  =  2y-3T,(y)  +  Ts(y), 

L2(y)  =  L5TJ{y)-2.5T!(y)  +  T7{y), 

where  T„(y)  is  the  Chebyshev  polynomial.22 

For  the  chosen  parameters,  the  matrices  A„,  B„,  and  C„  are  computed.  It  is  found  that, 
except  for  matrix  Au  all  other  matrices,  A„,  are  stable.  In  fact,  A\,  which  is  associated  to  wave 
number  do,  has  a  pair  of  unstable  eigenvalues  and  the  remaining  four  eigenvalues  of^i  are  stable. 
For  any  nonzero  initial  condition  the  solution  of  the  uncontrolled  system  diverges.  A  twelfth- 
order,  open-loop  system  including  the  dynamics  for  wave  numbers  a0  and  2a0  was  simulated  with 
the  initial  condition  «i(0)  =  6i(0)  =  1000(5-4,  4  05,  2-7)r,  and  a2( 0)  =  b2( 0)  =  ai(0)/0-18.  The 
divergent  and  oscillatory  responses  are  shown  in  figure  2.  The  plots  are  shown  for  the  velocity 
perturbations  ux(x* vx(x*y*,t)  at  the  point  (x*  y*)  =  (0.75^,-0.409),  and  also  for 

C/(x*,-l,/),  where  Cf  =  .5U^{diJdy)  is  evaluated  at  (x  =  x*  y*  =  0).  Although  a  control 

system  can  be  designed  to  modify  the  responses  of  a  set  of  modes,  for  simplicity,  stabilization  of 
dynamics  of  only  the  unstable  mode  is  considered  in  this  report.  The  spatial  distribution  of  the 
Lorentz  forces  is  assumed  as  given  in  equation  (21)  for  the  purpose  of  control. 

Simulation  results  are  obtained  using  an  optimal  control  system  and  a  simple  controller  that 
uses  wall-shear  stress  feedback  at  one  point.  For  simulation  of  the  optimal  controller,  it  is 
assumed  that  complete  information  on  the  state  variable,  xi,  of  the  subsystem  associated  with 
wave  number  a0  to  be  controlled,  is  available.  Unlike  the  optimal  controller,  synthesis  of  the 
simple  controller  is  relatively  easy,  because  only  a  measured  signal  at  the  wall  is  required  for 
feedback.  For  optimal  controller  synthesis,  each  component  of  vector  xx  must  be  known. 

5.1  OPTIMAL  CONTROL 

Using  the  optimal  control  technique,  a  feedback  controller  was  designed  for  the  decoupled 
system  in  equation  (22)  describing  the  dynamics  associated  with  the  fundamental  wave  number.  A 
quadratic  term  Oy]  was  used  in  the  performance  index  with  O  =  0.01  and  X  =  30.  The  weighting 
matrix  Q  and  X  were  chosen  after  several  trials  by  observing  the  simulated  responses.  The 
complete  closed-loop  system  including  the  dynamics  associated  with  wave  numbers  a0  and  2a0 
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and  the  feedback  control  law,  E  =  -Kxu  was  simulated.  The  initial  condition  of  figure  2  was 
retained.  It  was  assumed  that  longitudinal  control,  FUi  =  1;  surface-normal  control,  Fm  =  1;  and 
the  remaining  parameters  F\ik  in  equation  (21)  were  set  to  zero.  Selected  responses  are  shown  in 
figures  3  and  4  using  longitudinal  and  surface-normal  force  control,  respectively.  The  plots  are 
shown  for  the  electrode  voltage  E,  velocity  perturbations  u\{x*,y*,t),  V].(x* y*,t),  and  C/(x*, 0,t). 

In  the  closed-loop  system,  well-damped  responses  with  moderate  input  voltage  are  observed.  The 
corresponding  interaction  parameter  N*  =  /?(Lorentz  force/inertial  force)  is  0.0678 E,  where  the 
Lorentz  force  at  the  wall  is  45E(N)  and  the  inertial  force  is  (l l2)ju(dU/cty)  evaluated  at  the  wall. 
As  predicted,  the  surface-normal  force  control  requires  considerably  larger  electrode  voltage 
compared  to  the  longitudinal  force  control,  since  oto  <  -q.  It  is  important  to  note  that  using  this 
control  law,  stabilization  of  all  modes  is  accomplished,  even  ifB„,  n  *1  are  nonzero,  since  only 
exponentially  decaying  signal  X\  of  the  closed-loop  system  is  superimposed  on  the  stable  dynamics 
for  the  remaining  wave  numbers.  However,  estimation  of  signal  xx  from  measured  signals  must  be 
made  for  the  synthesis  of  the  optimal  controller. 

5.2  STABILIZATION  BY  WALL-SHEAR  STRESS  FEEDBACK 

Because  the  wall-shear  stress  is  a  constant  multiple  of  du^dy ,  a  simple  feedback  of  the 
derivative  of du\(x*y  =  0 ,t)ldy  is  considered  by  choosing  E  =  -Kdui/dy.  In  the  derivative  dujdy, 
the  contribution  of  only  the  fundamental  component  xj  of  the  state  vector  is  retained  and  the 
higher  harmonics  are  neglected.  The  gain  K  is  chosen  such  that,  in  the  closed-loop  system,  (. A\ - 
KB\C\)  is  Hurwitz.  Suitable  values  for  longitudinal  force  control  are  Fn2=  1  and  K  = 

0.05 (Ux/h) ;  and  suitable  values  for  surface-normal  force  control  are  Fm  =0  and  K  = 

0.05 [Ux/h)q.  The  feedback  gains  were  selected  after  several  trials  so  that  the  peaks  in  the 
electrode  voltage  are  not  too  large  during  the  transient  period.  The  feedback  gain  K  for 
stabilization  is  a  function  of  the  position  x*  where  the  wall-shear  stress  is  measured.  With  this 
choice  of  gains,  the  input  influence  matrix  B\  is  the  same  for  both  force  controls.  The  remaining 
parameters  Fuk  were  assumed  to  be  zero.  To  reduce  voltage  magnitude,  the  initial  conditions  were 
chosen  to  be  one-ninth  the  values  of  (a,(0),  b,(0')')  (/'  =  1,2)  assumed  in  section  5.1.  Selected 
responses  for  the  longitudinal  and  surface-normal  force  control  are  shown  in  figures  5  and  6. 

The  responses  for  vj,  and  C/  are  identical  for  the  longitudinal  and  surface-normal  force 
control.  Unlike  the  optimal  controller,  responses  are  poorly  damped  and,  in  spite  of  smaller 
perturbations  in  the  initial  condition,  considerably  large  electrode  voltage  is  required.  This  result 
should  not  be  a  surprise  because  the  optimal  controller  design  requires  complete  information  on 
Xi.  For  this  simple  control,  [dujdy)  feedback,  which  is  merely  a  linear  combination  of  state 
variables,  is  used.  Larger  peaks  in  the  input  voltage  are  observed  for  lateral  control  in  this  case 
also.  For  the  choice  of  spatial  distribution  of  Lorentz  forces  made,  stability  in  the  closed-loop 
system  is  preserved,  even  if  the  complete  measured  wall-shear  stress  including  higher  harmonics  is 
fed  back  because  the  dynamics  associated  to  all  wave  numbers  except  the  fundamental  one,  are 
stable  and  the  control  input  influence  matrix  Bn  =  0  for  n  =  0  and  n  >  1 .  These  responses  can  be 
improved  when  a  feedback  controller  that  uses  wall-shear  stress  signals  measured  at  several 
locations  on  the  wall  is  designed. 
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Figure  2.  Uncontrolled  System  at  x*  =  0. 75tj ,  y  *  =  -0.409  Above  the  Plate 
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Figure  6.  Surface-Normal  Force  Control  by  Wall-Shear  Stress  Feedback 

(x*  =  0. 751 ,  ,y*  =  -0.409) 
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6.  CONCLUSIONS 


This  report  examined  a  system-theory  approach  to  the  control  of  high  Reynolds  number  sea 
water  flow  over  a  flat  surface  using  electromagnetic  forces  produced  by  microtiles.  Navier-Stokes 
equations,  which  include  the  effect  of  Lorentz  forces,  were  linearized,  and,  using  the  Galerkin 
procedure,  a  linear,  state  variable  model  of  finite  dimension  was  derived.  Based  on  this  model, 
linear  feedback  control  laws  for  system  stabilization  were  obtained.  It  was  shown  that  by 
judiciously  selecting  spatially  distributed  Lorentz  forces,  asymptotic  stability  in  the  closed-loop 
system  can  be  accomplished  using  longitudinal  and  surface-normal  force  control.  For  lower  wave 
numbers,  longitudinal  force  control  is  more  effective  since  it  requires  smaller  electrode  voltage  to 
suppress  velocity  perturbations.  However,  a  control  system  that  uses  wall  surface  normal  force 
requires  smaller  electrode  voltage  compared  with  the  longitudinal  force  controller  when  control  of 
modes  associated  with  higher  wave  numbers  is  desired.  Although  stabilization  is  possible  using 
simple  feedback  of  measured  wall-shear  stress  signal,  an  optimal  controller  gives  faster  and  well- 
damped  responses.  However,  for  the  synthesis  of  the  optimal  controller,  it  is  necessary  to  design 
a  state  estimator  to  reconstruct  the  state  for  feedback  using  signals  measured  by  the  sensors. 

Several  problems  remain  to  be  solved  in  this  area:  questions  related  to  the  effect  of  mode 
truncation,  state  estimation,  control  and  observation  spillover,  digital  implementation,  and 
extension  of  the  system-theory  approach  to  nonlinearly  perturbed  dynamics,  which  plays  a  key 
role  in  turbulent  flow,  are  all  important  and  require  further  investigation. 
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